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Abstract. It was shown by Faltings jFal84| and Hriljac |Hri85j 
that the Neron-Tate height of a point on the Jacobian of a curve 
can be expressed as the self-intersection of a corresponding divisor 
on a regular model of the curve. We make this explicit and use it to 
give an algorithm for computing Neron-Tate heights on Jacobians 
of hyperelliptic curves. To demonstrate the practicality of our 
algorithm, we illustrate it by computing Neron-Tate heights on 
Jacobians of hyperelliptic curves of genus 1 < 5 < 9. 



1. Introduction 

The problem considered in this paper is that of computing the Neron- 
Tate (or canonical) height of a point on the Jacobian of a curve of 
genus greater than 2. For curves of genus 1 and 2 the existing methods 
(classical in genus 1, and due to Flynn, Smart, Cassels and others 
in genus 2 - see for example |CF96j and |FS97j ) make use of explicit 
equations for projective embeddings of Jacobians, and have proven to 
be very successful in practice. It does not seem practical at present 
to give explicit equations for Jacobians of curves of genus 3 and above 
(see |Stu01j and |MuelOj for recent attempts and an examination of the 
difficuties faced). We propose an alternative approach to computing 
the Neron-Tate height based on Arakelov theory. To demonstrate that 
our method is practical, we give numerical examples where we compute 
heights of points on Jacobians of hyperelliptic curves of genus I < g < 
9. 

The main application of these computations is at present the compu- 
tation of regulators of hyperelliptic Jacobians up to rational squares; 
this will allow the verification up to rational squares of the conjec- 
tures of Birch and Swinnerton-Dyer. The author has also developed an 
algorithm to bound the difference between the naive and Neron-Tate 
heights, again using Arakelov theory [HollOj . Together, the algorithms 
allow a number of further applications such as computing a basis of 
the Mordell-Weil group of a hyperelliptic Jacobian, computing integral 
points on hyperelliptic curves (see |BMS"'"10] ) and of course verifying 
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the conjectures of Birch and Swinnerton-Dyer, up to the order of the 
Shafarevich-Tate group. 

Whilst the theoretical sections of this paper are largely indepen- 
dent of the curve chosen, the very geometric nature of Arakelov theory 
means that the details of the algorithm, and especially its implemen- 
tation, will depend greatly on the geometry of the curve considered. 
Let C be a curve defined over a number field k. Our method makes a 
number of computational assumptions: 

(a) We have a uniform and convenient way of representing divisor 
classes on the curve C. 

(b) We are able to rigorously compute abelian integrals on C to 
any required (reasonable) precision. 

(c) We are able to write down a regular proper model for C over 
the integers Ok (though this need not be minimal). 

(d) For each non-Archimedean place u that is a prime of bad re- 
duction for we are able to compute the intersection matrix 
of the special fibre 

(e) We have a way of computing Riemann-Roch spaces of divisors 
on C. 

Assumptions (a) and (b) cause us to restrict our attention to hyper- 
elliptic curves. The computer algebra package MAGMA |BCP97] has in- 
built commands to deal with all of these for hyperelliptic curves over 
number fields. For (a), it uses Mumford's representation (see |MM84t 
3.19]) for divisor classes. A MAGMA implementation by P. van Wamelen 
is used for integral computations in (b); this does not use the usual 
numerical integration techniques as these are inherently non-rigorous; 
instead, hyperelliptic functions are locally approximated by truncated 
power series and formally integrated. The computation of the inter- 
section matrices of the special fibres at the bad places is produced by 
MAGMA as by-product of the computation of the regular proper model, 
implemented by S. Donnelly using techniques as in |Liu02| Chapter 8]. 
For computing Riemann-Roch spaces, MAGMA makes use of the method 
of Hess |Hesj . 

In order to simplify the exposition, we restrict our attention to curves 
with a rational Weierstrass point. As such, unless otherwise stated, C 
will denote an odd-degree hyperelliptic curve over a number field k, 
and ^ a proper regular — though not necessarily minimal — model of C 
over the integers Ok- 

The author wishes to thank Samir Siksek for introducing him to this 
fascinating problem, as well as for much helpful advice and a careful 
reading of this manuscript. Thanks are also due to Martin Bright and 
Jan Steffen Miiller amongst others for very helpful discussions. 
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2. A Formula of Faltings and Hriljac 

As before C is an odd degree hyperelliptic curve over a number field 
k. We fix once and for all the following notation: 

• Mk a proper set of places of k\ 

• the subset of non-Archimedean places of k; 

• the subset of Archimedean places of k; 

• k(z/) the residue field at a z/ G M°; 

• iy the usual intersection pairing between divisors over u G M° 
(see |Lan88l IV, §1]). 

We shall make use of the following result which can be found in 
Lang's book |Lan88[ IV, §2]. 

Theorem 1. (Faltings and Hriljac) Let D be a degree zero divisor on 
C , and let E he any divisor linearly equivalent to D hut with disjoint 
support. Then the height with respect to the -d- divisor of the point on 
Jac(C) corresponding to D is given hy 

where $ and gD,u are defined as follows: 

- $ sends a divisor on the curve C to an element of the group of 
fihral Q-divisors on ^ with order zero along the irreducihle component 
containing infinity, such that for any divisor D on C (with Zariski 
closure D ) and fihral divisor Y on 'is , we have Lu {^D + $(-D), F) = 0. 

- gD^i, denotes a Green function for the divisor D, when C is viewed 
as a complex manifold via the emhedding v. 

Our strategy for computing the Neron-Tate height of a degree zero 
divisor D using the above formula is as follows: 

(1) Determine a suitable divisor E as above. This is explained in 
Section [31 

(2) Determine a finite set £H C M° such that for non-Archimedean 
places V not in we have ly (^D + ^{D), E) = 0. This is 
explained in Section HJ 

(3) Determine l,^ for z/ G This is explained in Sec- 
tion o 

(4) Determine L,y (^D, !E) for u E Ul. This is explained in Section El 

(5) Compute the Green function g£) y{E) for Archimedean v. This 
is explained in Section [71 

In the final section we give a number of worked examples. 

By a straightforward Riemann-Roch computation, we can write down 
a divisor in Mumford form that is linearly equivalent to D. We replace 
D by this Mumford divisor. Thus we may suppose that D = D' — d-oo 
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where d < g and D' = zeTos{a{x),y — b{x)) with a(x), b{x) G k[x] 
satisfying certain conditions as given in jMM84t 3.19, Proposition 1.2]. 

3. Step 1: Choosing E 

If the support of D' does not contain Weierstrass points, choose A G /c 
such that a(A) 7^ 0, and set 

(1) E = 'mv{D') — -zeros(a; — A) , 

where inv denotes the hyperelhptic involution. If d is odd, this is not 
a divisor but a Q-divisor. This is unimportant, but if the reader is 
troubled he or she should multiply E hj 2, and then appeal to the 
quadraticity of the height. 

If the support D' does contain Weierstrass points, either: 

a) replace D by a positive multiple of itself to avoid this, or 

b) add a divisor of order 2 to to remove them. 

(a) is simpler to implement, (b) generally faster computationally. 
Now E ~iin — -D, and so 

h^{0{D))= J2^os\<i^)\'u(D + ^D),E)+^ Yl 9dAE). 

This is seen by viewing the expression on the right hand side as the 
global Neron pairing on divisor classes, which is a quadratic form; since 
E is linearly equivalent to —D a minus sign results, which cancels with 
those in Theorem [1] to yield the above expression. 

4. Step 2: Determining a Suitable 9^ 
We wish to find a finite set £H C M° such that 

(2) 6, + =0 

for all z/ ^ £H. To make this as general as possible, we will for the mo- 
ment just assume C is a smooth curve over k in the weighted projective 
space Pfc(ao! • • • ? o.n)- Let C denote its closure in Pc'^(ao, . . . , ctn). Let 
Qi denote the set of places of bad reduction for C", outside which C 
is smooth over Ok- 

It suffices to solve our problem for prime divisors, as we can then ob- 
tain results for general D and E easily. Let X and Y be prime divisors 
on C , and let d be the degree of Y . Let Hi, ... , Hd+i be a collection 
of weighted integral homogeneous forms of degrees ei,...,ed+i > 
on Poj.(ao, . . . , ci„), geometrically integral on the generic fibre and co- 
prime, such that for all pairs i 7^ j, we have on the generic fibre that 
Hi n Hj n C = $ (we will confuse Hi with the hypersurface it defines) . 

Let Q2 be the set of z/ G M° \ Qi such that 

(H,), n (Hj), n c: ^ 0. 
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for some I < i, j < d + 1. 

Let FF{Y) denote the function field of Y, a finite extension of k, 
and let Ny : FF{Y) — )■ k denote the norm map. In practice, we can 
use Grobner bases to find an isomorphism FF{Y) — )■ k[t]/a{t), and so 
can readily compute Ny. 

Let be integral weighted homogeneous equations for X 

such that no fi vanishes on Y, and set deg(/j) = dj. Finally let Q3 be 
the set of z/ G M° \ (Qi U Q2) such that 

ord, (Ny (/JVi/f^)) ^0 

for some l<j<r, l<'j<(i + l. 

Lemma 2. Set 

'K = QiUQ2UQ3. 
Ifvi'K then (X + $(X),F) =0. 

Proof. Outside Qi, C is smooth over C^, and hence it is regular and 
all its fibres are geometrically integral. As a result, 

(3) (X + F) = (X, Y) for z/ ^ ^. 

Suppose (X,F) ^ 0, so (X)^ n (F)^ ^ 0. We will show z/ e 9=1. 

Recall that X and Y are cycles on C of relative dimension zero over 
Ok, and so their fibres over closed points are cycles of dimension zero. 
Observe that, since the fj are integral, we must have zeros (/j) D X for 
all j . Hence there is some jo such that /jp vanishes on some irreducible 
component of (equivalently, closed point in) {Y)^, (in fact, this holds 
for any jo)- As a result, Lu (zerosc/(/jo) , F) > 0, since we assume 
zerosc"(/jo) and Y have disjoint support on the generic fibre. Now 
suppose 1/ ^ Then for all 2, since u ^ Q3, we must have 

ord. (Ny (/l;//^^))=0. 

Hence by [LanSSj HI, Lemma 2.4, p56], we see that for all i, 

o-.(div.(|,),F) 

= Ci- (zerosc7'(/jo) > ^) - c?io " (zerosc'(i^i) , Y) . 
Now as Cj > and dj^ > 0, we see that for all i, 

(4) [zeiosciHi) ,Y) > 0, 

so every zerosc"(-f^i) meets Y. But since v ^ Q2 the zero-dimensional 
cycles zerosc"(-ffi) HF are pairwise disjoint. Hence the zerosc"(-ffj) ny 
cannot all meet the zero- dimensional cycle {Y)y as it has degree d\ 
deg(F) = deg{{Y)y) = d. This contradicts Equation HI and so we are 
done. □ 
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Note that if Oq = . . . = a„ = 1, so C hves in 'unweighted' projec- 
tive space, then the Hi should all be taken to be hyperplanes and the 
argument simplifies somewhat. 

5. Step 3: Determining ij. ($(D), £') 

We next discuss the computation of the term for a non- 

Archimedean place V. Recall that by our assumptions in Section |2] we 
are able to write down a proper regular model and the intersection 
matrix for 'lou for all bad places u. Clearly vanishes if 

is integral, in particular if z/ is a good prime. Suppose i/ is a bad prime. 
Since we have the intersection matrix of "iou-, the problem is reduces to 
the following: 

Problem 3. Given a finite place v, a horizontal divisor X and a prime 
fibral divisor Y over u, compute l^, {X, Y) . 

We may replace the base space S = Spec(Cfc) by its completion S at 
u ( [Lan88l III, Proposition 4.4, page 65]), and we may further assume 
that X is a prime horizontal divisor on ^ x 5 S*. By Lemma S] below, 
this means that the support of X^ is a closed point of and so we 
can find an affine open neighbourhood U = Spec (A) of X^, in ^ x 5 S*. 

Lemma 4. If X is a prime horizontal divisor on an arithmetic surface 
3^ over a p-adic local ring [O^v), then the support of Xy is a prime 
divisor on (in other words, X^, is irreducible hut not necessarily 
reduced). 

Proof. There exists a number field L/k and an order R in L such that 
X is isomorphic to Spec(i?). Write L = k[t]/a(t), where a monic and 
irreducible with integral coefficients. Let k, denote the residue field of 
k, and a the image of a in K[t]. If X^, is not irreducible, then there 
exist f,g& K[t] coprime monic polynomials such that f ■ g = a. This 
factorization of a lifts to a factorization of a by Hensel's Lemma. □ 

Now it is easy to check whether X^ fl F = 0; if so, l^, (X, Y) = 0. 
Further, if Xi, C Y and X^, is not contained in any other fibral prime 
divisor, then tj, {X,Y) = deg(X); this is easily seen since locally Y = 
zeroSi7(z/), and we can take the norm of u from the field of fractions 
FF{X) down to k. 

We are left with the case where X^ lies at the intersection of several 
fibral prime divisors. We find equations fi,...,f^EA for Y 

as a subscheme of U^. Then choose any lifts fi of /j to A. Now we 
need two easy results in commutative algebra: 

Lemma 5. let R he a ring, p (z R any element, and I an ideal contain- 
ing p. Suppose we have ti, . . . ,tr & R such that the images ti, . . . ,tr in 
R/{p) generate the image of I inR/{p). Then I = {ti, . . . ,tr,p) . 
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Proof. Let x G /. Write x for the image of x in R/{p), and write 
X = ^i^ittA for some a G R/{p)- Choose hfts of a, to i?. Then 

y =^ X — ^[=1 ajtj hes in p ■ i?. Hence x is in (tj, . . . ,tr,p) and so 
/ C {ti, . . . ,tr,p). Now p G / by assumption, and ti G I/{p), so there 
exists (yfj in J such that Qi — ti ^ p ■ R, so & I . □ 

Lemma 6. Let R be a regular local ring, and ti, . . . ,tr G R be such 

def 

that I = (ti, . . . ,tr) is a prime ideal of height 1. Now I is principal; 
write I = (t). Then there exists an index i and a unit u E R such that 
ti = tu. In particular, there exists an index i with I = (ti). 

Proof. Ris a unique factorization domain, and so for each i we can write 
ti = t'^t for some t- G R. Hence I = t ■ {t[, . . . ,t'^), and {t[, . . . ,t'^) = 1. 
We want to show some t- is a unit. Suppose not; then since A is local, 
all the t'i lie in the maximal ideal, so {t[, . . . ,t'^) is contained in the 
maximal ideal, a contradiction. □ 

Now from these we see that one of the fi or u must be an equation 
for y in a neighbourhood of (and it cannot be as Xi, lies on an 
intersection of fibral primes). Now if any fi vanishes on X^,, it cannot 
be the fi we seek. Exclude such fi, and then for each of the remaining fi 
compute its norm from FF{X) to the completion of k. The minimum 
of the valuations of such norms will be achieved by any fi which is a 
local equation for Y at X^,, and hence {Y, X) is equal to the minimum 
of the valuations of the norms. 

6. Step 4: Determining l,, {D,'E) 

Finally, we come to what appears to be the meat of the problem 
for non-Archimedean places: given two horizontal divisors D and E 
and a place z/ G fH, compute the intersection (D, 'E) . However, the 
techniques used in previous sections actually make this very simple. 

Fix a non- Archimedean place v. Let 5* denote the z/-adic completion 
of 5*, and set ^ = ^ S'. It is sufficient to compute the intersection 
Ly (X, Y) where X and Y are prime horizontal divisors on in par- 
ticular (by Lemma H]), the supports of Xy and Y^ are closed points of 

Now if Supp(X,^) 7^ Supp(F,^), then Ly{X,Y) = 0. Otherwise, let 
U = Spec(A) be an affine open neighbourhood of Supp(Xy). Let 
fi, . . . , fr generate the ideal of X on U; then by Lemma [6] we know 
that some fi generates the ideal of X in a neighbourhood of X,y. If fj 
vanishes on Y, we can throw it away. We obtain 

Proposition 7. 

(X,y) =min (ord,(/ar])) 

i 

as i runs over {1, . . . , r} such that fi does not vanish identically on Yy. 
Here f\Y] is defined to be either 
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(1) the norm from FF{Y) to the completion of k of the image of fi 
in FF(Y), or equivalently, 

(2) Ylj fiiPj)^^ where Y = '^jTijPj over some finite extension l/k 
(see [Lan88l II,§2, page 57];. 

Proof. If fi is not identically zero on Y"^, then zeros<^-;(/i) and Y have no 
common component and moreover fi is regular on a neighbourhood of 
Y so Li, (poles<:^(/i), F) = 0, and so [LanSSj II, Lemma 2.4, p56] shows 
that 

(5) (zeros.^</,) , Y) = ord, (/.[F]) . 

Now zeroSc^(/j) > X, so ord^(/j[F]) > Li,{X,Y). Moreover, by 
lemma O there is an index i^ such that /jg generates X near Xj^, and 
since = Y^ is irreducible we have that 

iv (X, Y) = (zeros^-(/i J , F) = V lengthy ( , f ^ ) 

(6) ^ V(/.) + /W 

= length^^^^ 



O 



Y 

where the sum is over closed points p oi ^ lying over u, and Jy is the 
defining ideal for Y in the local ring under consideration. Now any 
other fi will have 

(7) (zeros^~(/i) ,Y) > (zeroSc^(/iJ , Y) , 

so the result follows. □ 

As regards the computation of the fi[Y], definitions (1) and (2) given 
in Proposition [7| lead to slightly different approaches, but both make 
use of Pauli's algorithms |PR01j . In our implementation, discussed in 
Section [HI we use (2) as it seems easier; however (1) may lead to an 
implementation that is faster in practice. 

7. Step 5: Computing goA^) 

Finally, we must compute the Archimedean contribution. Fix for the 
remainder of this section an embedding a of in C corresponding to a 
place u G M^. Let C^- denote the Riemann surface corresponding to 

7.1. The PDE to be Solved. As a starting point, we take |Lan83l 
Chapter 13, Theorem 7.2], which we summarise here. 

Given a divisor a on of degree zero, let a; be a differential form on 
Ca such that the residue divisor res(a;) equals a (such an u can always 
be found using the Riemann- Roch Theorem). Normalise u by adding 
on holomorphic forms until the periods of u are purely imaginary. Let 

(8) dga u + u. 
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Then ga is a Green function for a. Thus it remains to find, normahse 
and integrate such a form u. 

7.2. Application of theta functions to the function theory of 
hyperelliptic curves. We can use "i^-functions to solve the par- 
tial differential equation (|8]) of Section I7.1[ in a very simple way. For 
background on ?9-functions we refer to the first two books of the 'Tata 
lectures on theta' trilogy, |Mum83] . |MM84] . '(^-functions are complex 
analytic functions on which satisfy some quasi-periodicity condi- 
tions, thus they are an excellent source of differential forms on the 
(analytic) Jacobian of C^. To get from this a differential form on C^r 
we simply use that C„ is canonically embedded in Jac(Co-) by the Abel- 
Jacobi map, so we can pull back forms from Jac(Co-) to C„. 

Fix a symplectic homology basis Ai, Bi on C^^ as in |MM84j : by this 
we mean that if i(— , — ) denotes the intersection of paths, then we 
require that the Ai, Bi form a basis of Hi{Ccr, Z) such that 

i{Ai, Aj) = i{Bi, Bj) = for z ^ j 

and 

i{Ai,Bj) = Sij. 

We also choose a basis ui, . . .Ug of holomorphic 1-forms on C^, nor- 
malised such that 




We recall the definition and some basic properties of the multivariate 
'i^-function: 

(9) -^[z] Q) =^ exp{ninQn^ + 27rm ■ z) 

nin Zs 

which converges for z in and Q a g x g symmetric complex ma- 
trix with positive definite imaginary part. The function satisfies the 
following periodicity conditions for rn,n in Z^: 

(10) ^{z + m;n) =^{z;n), 

(11) '&{z + nfl; il) = exp(—ninp,n^ — 2'Kinz) d{z] Vl). 

We will set Vt to be the period matrix of the analytic Jacobian of Ca 
with respect to the fixed symplectic homology basis (as in |MM84] ) . 
and z will be a coordinate on the analytic Jacobian. This means that 

'i 

Let 

/II 1 1\ 1 . 
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V2' 2 ' ' '27 2 
A = n-S' + 5". 

Then |MM84l Theorem 5.3, part 1] tells us that ^9(A - z) = if and 
only if there are Pi, . . . Pg-i in such that 

z 



= J] / w (mod + 1]Z»). 



This is a crucial result which allows us to construct a quasifunction on 
Jac(Co-) with prescribed zeros, and from this obtain the Green function 
we seek. 

7.3. Solution of the Partial Differential Equation. Let D, Dq be 

two effective divisors of degree g on with disjoint support, containing 
no points at infinity, nor any pairs p + q oi points such that p = inv(g). 
Then the classes [0{D — g ■ oo)] and [0{Dq — g ■ oo)] lie outside the {}- 
divisor on the Jacobian; indeed, the association D i— > [0{D — g ■ oo)] is 
an isomorphism from divisors with the above properties to Jac(Co-) \ 
see |MM84l 3.31]. Write a : Div(C^) Jac(C^) for the map sending 
a divisor E to the class [0{E — deg{E) ■ oo)]. 
For z in Jac(Co-) we set 



^{z + A-a{Do)y 
Then for p in we set E{p) = G{a{p)) so 
(12) F(p) ''("W + A-„(fl)) 



tf(a(p) + A-a(D„))- 



If we let CO = dlogF{p) then it is clear that res(w) = D — Dq. It 
then remains to normalise u to make its periods purely imaginary, and 
then integrate it. We have a homology basis Ai, Bi, and we find: 

uj= d\ogF{p) = \ogG{a{p) + e^) - \ogG{a{p)) = 

JAk JAk 

(where Cfc = (0, 0, . . . 0, 1, . . . , 0) with the 1 being in the k-th. position), 
and 

uj= dlogF{p) =logG{a{p) + n.ek) -\ogG{a{p)) 

= 27riel -{aiD) -a{Do)). 
From this we can deduce that the normalisation is 
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uj — dlog 



i9{a{p) + A - a{D)) 



- 27ri [(Im(l]))-^ Im(a(D) - a{Do))] 



0J2 



00, 



9 J 



where p is a point on C„. 



Now we integrate to get the Green function gn-Doip) = /<^^ O' + O', 



where ooc„ denotes the point at infinity on €„: 



QD-Doip) = 2 log 



+ 47r [(Im((]))-i Im(a(D) - a(Do))] . Im 



V 



0J2 



UJg jj 



^ ^{a{p) + A - a{D)) 

^ ^{a{p) + A-a{D,)) 
+ 47r(Im(fi))-^ • lm{a{D) - a{Do)) ■ Im {a{p)) . 

Given divisors Dq and Eq containing no Weierstrass points or 
infinite points or pairs of points which are involutions of each other, and 
having disjoint support, we can use this formula to compute \gD-Do\E— 
Eq\ which is simply the product over points p G Supp(£^ — £^0) of the 
complex absolute value of g{p). We are done. 

8. Examples 

We have created a test implementation of the above algorithm in 
MAGMA. The following results were obtained using a 2.50 GHz Intel 
Core2 Quad CPU Q9300: 

First, we let C/Q be the genus 3 hyperelliptic curve given by 

C -.y"^ = x^ - 15x^ + llx^ - 13a; + 25. 

Let E he the points on the Jacobian corresponding to the degree 
divisors (1, 3) — 00, = (0, —5) — 00 respectively. We obtain the 
following: 

h{D) = 1.77668... 
h{E) = 1.94307... 
h(D + E)^ 4.35844 . . . 
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h{D-E) = 3.08107... 

2h{D) + 2h{E) - h{D + E) - h{D - E) = 1.26217 x lO'^^ 

with a total running time of 31.75 seconds. We note that our result is 
consistent with the parallelogram law for the Neron-Tate height, which 
provides a useful check that our implementation is running correctly. 

Next we give two families of curves of increasing genus. Firstly the 
family y"^ = x^^^^ + 2x'^ — lOx + 11 with D denoting the point on the 
Jacobian corresponding to the degree divisor (1, 2) — oo (all times are 
in seconds unless otherwise stated): 
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h{D) 


time 


1 


1.11466. . . 


1.94 


2 


1.35816... 


6.44 


3 


1.50616... 


15.10 


4 


1.61569... 


32.71 


5 


63.4292 . . . 


72.23 


6 


1.77778... 


212.37 


7 


51.0115... 


20 minutes 


8 


1.89845 . . . 


3 hours 


9 


78.8561 . . . 


16 hours 



Now we consider the family y"^ = x'^^~^^ + 6x'^ — 4x + l with D denoting 
the point (1, 2) + (0, 1) — 2 ■ oo on the Jacobian: 
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time / seconds 


1 


1.41617. . . 


2.06 


2 


1.37403 . . . 


6.73 


3 


1.50396 . . . 


15.62 


4 


1.40959 . . . 


32.60 


5 


1.70191... 


76.48 


6 


1.81093... 


291.17 


7 


1.71980... 


1621.50 



A fully-functioning and more efficient implementation of Neron-Tate 
height computations is currently being carried out in MAGMA by J. S. 
Miiller, combining ideas from this paper with those from his own PhD 
thesis, [MuelOj . Miiller's approach to computing Li,{D,E) is quite 
different from that used here; he uses Grobner bases to compute directly 
the Op-length of the modules j^^j^ as p runs over closed points of the 
special fibre, in contrast to the method in this paper where we compute 
norms down to the ground field and then compute valuations. 
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